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Abstract We investigate observational constraints on the Dvali, Gabadadze and Porrati 
(DGP) model with Gamma-ray bursts (GRBs) at high redshift obtained directly from the 
Union2 Type la supernovae data (SNe la) set. With the cosmology-independent GRBs, 
the Union2 set, as well as the cosmic microwave background (CMB) observations from 
the WMAP7 result, the baryon acoustic oscillation, the baryon mass fraction in clusters 
and the observed H(z) data, we obtain that the best-fit values of the DGP model are 
{fl M0 , n r J = {0.235^;^|, 0.138±g;^}, which favor a flat universe; and the transi- 
tion redshift of the DGP model is zt = 0.67^ 04,- These results lead to more stringent 
constraints than the previous results for the DGP model. 

Key words: Gamma rays : bursts — Cosmology : cosmological parameters 
1 INTRODUCTION 

The accelerating expansion of the current universe has been confirmed by recent cosmological observa- 
tions, such as Type la supernovae (SNe la; Riess et al. 1998; Perlmutter et al. 1999; Amanullah et al. 
2010), cosmic microwave background (CMB; Bennett et al. 2003; Spergel et al. 2003; Komatsu et al. 
2010), large scale structures (LSS, Tegmark et al. 2004; Eisenstein et al. 2005), as well as the x-ray gas 
mass fraction of clusters (Allen et al. 2004). By assuming General Relativity, a dark energy component 
with negative pressure in the universe has been invoked as the most feasible mechanism for the accel- 
eration. In addition to the cosmological constant (the ACDM model), many candidates of dark energy 
have been taken into account. Examples include the scalar field models with dynamical equation of 
state [e.g., quintessence (Ratra and Peebles 1988; Caldwell et al. 1998; Choudhury and Padmanabhan 
2005), phantom (Caldwell 2002; Wu and Yu 2006), quintom (Feng et al. 2005; Guo et al. 2005; Liang et 
al. 2009), k-essence (Armendariz-Picon et al. 2001; Chiba 2002), tachyon (Padmanabhan 2002; Frolov 
et al. 2002)], the Chaplygin gas (Kamenshchik et al. 2001) and the generalized Chaplygin gas model 
(GCG, Bento et al. 2002; Zhu 2004), the holographic dark energy (Cohen 1999; Li 2004), the agegraphic 
dark energy (Cai 2007; Wei & Cai 2008), the Ricci dark energy (Gao et al. 2009) and so on. 

On the other hand, many alternatives to dark energy in which gravity is modified have been pro- 
posed as a possible explanation for the acceleration. Examples include the f(R) theory in which the 
Einstein-Hilbert action has been modified (Capozziello & Fang 2002; Vollick 2003; Carroll et al. 2004); 
the Cardassian expansion model in which the Friedmann equation is modified by adding an extral 
Cardassian term (Freese and Lewis 2002; Wang et al. 2003; Zhu and Fujimoto 2002, 2003); as well 
as the braneworld models, in which our observable universe is considered as a brane embedded in a 
higher dimensional bulk spacetime and the leakage of gravity force propagating into the bulk can lead 
to the current accelerated expansion of the universe (Randall and Sundrum 1999). 
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In 2000, Dvali, Gabadadze and Porrati proposed a 5-dimensional brane world model in which a self- 
accelerating branch is included (the so-called DGP model, Dvali et al. 2000). The dynamics of gravity 
is governed by a competition between a Ricci scalar term in the 4-dimensional brane and an Einstein- 
Hilbert action in the 5-dimensional bulk. The Friedmann equation of the DGP model is modified as 



where H is the Hubble parameter as a function of redshift z, SImo an d represent the fractional 
contribution of the matter and curvature, and fl Ic = I/At^Hq is the bulk-induced term respect to 
the crossover radius r c . For scales below r c , the induced 4-dimensional Ricci scalar dominates and the 
gravitational force is the usual 1/r 2 behavior; whereas for distance scales larger than r c , the gravitational 
force follows the 5-dimensional 1/r 3 behavior. The normalization condition can be given by £Ik + 

(yfi^ + y/£l Ic + f^Mo) = 1; for a spatially flat scenario, tt Ic = (1 — S1mo) 2 /4. 

The DGP model is a testable scenario with the same number parameters as the standard ACDM 
model and has been constrained from many observational data, such as SNe la (Deffayet et al. 2002; 
Avelino and Martins 2002; Zhu and Alcaniz 2005; Maartens and Majerotto 2006; Barger et al. 2007; 
Reboucas 2008), the angular size of compact radio sources (Alcaniz 2002), the baryon mass fraction 
in clusters of galaxies (CBF) from the x-ray gas observation (Zhu and Alcaniz 2005; Alcaniz and 
Zhu 2005), CMB (Lazkoz et al. 2006; Rydbeck et al. 2007; He et al. 2007), the large scale structures 
(Multamaki et al. 2003; Lue et al. 2004; Koyama and Maartens 2006; Song et al. 2007) and the baryon 
acoustic oscillation (BAO) peak (Guo et al. 2006), the observed Hubble parameter H (z) data (Wan et al. 
2007), the gravitational lensing surveys (Jain et al. 2002; Zhu and Sereno 2008), the age measurements 
of high-z objects (Alcaniz, Jain and Dev 2002) and the lookback time to galaxy clusters (Pires, Zhu and 
Alcaniz 2006); as well as some different combined data (Bento et al. 2006; Davis et al. 2007; Movahed 
et al. 2009; Xia 2009; Li et al. 2010). See Lue (2006) for review on the DGP phenomenology. 

Recently, Gamma-ray bursts (GRBs) have been proposed as distance indicators and regarded as a 
complementary cosmological probe to the universe at high redshift (Schaefer 2003; Bloom et al. 2003; 
Dai et al. 2004; Ghirlanda et al. 2004; Friedman and Bloom 2005; Firmani et al. 2005, 2006; Liang 
and Zhang 2005; Bertolami and Silva 2006; Ghirlanda et al. 2006; Schaefer 2007; Wright 2007; Wang 
et al. 2007; Amati et al. 2008; Basilakos and Perivolaropoulos 2008; Mosquera Cuesta et al. 2008a, 
2008b; Qi et al. 2008a, 2008b). For constraints on the DGP model from GRBs with their associated 
joint observations, see some recent works, e.g., Wang et al. (2009a); Wei (2010a); Xu and Wang (2010). 
However, the empirical luminosity relations of GRBs have usually been calibrated by assuming a cer- 
tain cosmological model with particular model parameters, due to the lack of the low-redshift sample. 
Therefore the calibration are always cosmology-dependent and the so-called circularity problem occurs 
in GRB cosmology. The circularity problem cannot be avoided completely by means of statistical ap- 
proaches (Schaefer 2003; Li et al. 2008; Wang 2008; Samushia and Ratra 2010; Xu 2010), because an 
input cosmological model is still required. Liang et al. (2008) presented a new method to calibrate GRB 
luminosity relations in a completely cosmology-independent way: GRB sample in the redshift range of 
SNe la are enough to calibrate GRB relations and their luminosity distances can be obtained directly 
from SNe la by the interpolation method or by other similar approach (Liang and Zhang 2008; Kodama 
et al. 2008; Cardone et al. 2009; Gao et al. 2010; Capozziello and Izzo 2010). Following the cosmology- 
independent GRB calibration method, the derived GRB data at high redshift can be used to constrain 
cosmological models by using the standard Hubble diagram method (Capozziello and Izzo 2008; Izzo 
et al. 2009; Wei and Zhang 2009; Wei 2009; Qi et al. 2009; Wang et al. 2009a, 2009b; Wang and Liang 
2010; Liang, Wu and Zhang 2010; Liang, Wu and Zhu 2010; Wei 2010a, 2010b; Freitasa et al. 2010; 
Liang, Xu and Zhu 2010; Demianski et al. 2010). 

Very recently, Liang, Wu and Zhu (2010) calibrated GRB data at high redshift directly from the 
Union2 compilation of 557 SNe la data set (Amanullah et al. 2010); and constrained the Cardassian 
model (Liang, Wu and Zhu 2010) and the generalized Chaplygin gas (GCG) model (Liang, Xu and Zhu 
2010) by combining the updated GRB data with the joint observations, such as the Union2 set of SNe la, 
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the CMB observation from the seven-year data of Wilkinson Microwave Anisotropy Probe (WMAP7; 
Komatsu et al. 2010) result and the BAO observation from the spectroscopic Sloan Digital Sky Survey 
(SDSS) galaxy sample (Eisenstein et al. 2005). 

In this paper, we investigate observational constraints on the DGP model including the updated the 
distance moduli of the GRBs at high redshift obtained directly from the Union2 set. We combine the 
GRB data with the joint observations such as the Union2 set, the CMB observation from the WMAP7 
result; the BAO observation from the spectroscopic SDSS galaxy sample (Eisenstein et al. 2005); the 
baryon mass fraction in clusters of galaxies from the x-ray gas observation (Allen et al. 2004); and the 
observed Hubble parameter data (H(z); Simon et al. 2005; Gaztanaga et al. 2009). We also obtain the 
transition redshift zt of the DGP model. We find that the combination of these recent data sets tighter 
constraints on the DGP model, which favors a flat universe. The paper is organized as follows. In section 
2, we introduce the analysis for the observational data including the updated cosmology-independent 
GRBs, as well as the Union2 SNe la set, the CMB observations from the WMAP7 result, in addition to 
the BAO, CBF and H (z) data. In section 3, we present results which put constraints on the DGP model 
from the joint observations. Conclusions and discussions are given in section 4. 

2 OBSERVATIONAL DATA ANALYSIS 

The recent Union2 compilation consists of a 557 SNe la data set (Amanullah et al. 2010). In this paper, 
we use the updated the distance moduli of the 42 GRBs at z > 1.4 (Liang, Wu and Zhu 2010), which 
are obtained by the five luminosity relations (Schaefer 2007) calibrated with the sample at z < 1.4 by 
using the linear interpolation method from the Union2 set. For more details about the calculation, see 
(Liang et al. 2008; Liang, Wu and Zhang 2010). Constraints from SNe la and GRB data can be obtained 
by fitting the distance moduli fx(z). A distance modulus can be calculated as 

/' = 51o 8'r7^ + 25 = 51 °gio- D L "Mo, (2) 
Mpc 

where fj,o = 5 log 10 [7?o/(100km/s/Mpc)] + 42.38, and the luminosity distance Dl is calculated by 

dz' 



D L = H d L = (1 + z)Q k 1/2 sinn 



n 1/2 



E(z' 



(3) 



where sinn(a;) is sinh for fik > 0, sin for < 0, and x for Cl^ = 0. The \ 2 value of the observed 
distance moduli can be calculated by 

2 _ [MobsQi) ~ Kzj)} 2 

i=l M,' 

where fJ. bs(zi) are the observed distance modulus for the SNe la and/or GRBs at redshift Zj with its 
error er^; ^t(zi) are the theoretical value of distance modulus from cosmological models. Following an 
effective approach (Nesseris and Perivolaropoulos 2005), we marginalize the nuisance parameter fio by 
minimizing 

xl = C - B 2 /A, (5) 

where ,4 = E 1 /^.' 5 = £ [MobsOi) -51og 10 Ac]/^., and C = £ \}i ohs {zi) -51og 10 L> L ] 2 /cr^. 

For the CMB observation from the WMAP7 result (Komatsu et al. 2010), the shift parameter is 
constrained to be R = 1.725 ± 0.018, which can be expressed as (Bond et al. 1997) 



it = i 2j^ i l k ' sinn 



"k 



dz 



o 



(6) 



E(z 

is the redshift of recombination which is given by (Hu and Sugiyama 1996) 

1048[1 + 0.00124(fVi 2 r - 738 (l + gi(n M0 h 2 ) 92 )}, (7) 
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XCMB - n mo2 ' 



where g x = 0.0783(fM 2 )~ a238 (l + 39.5(0 6 /i 2 )-°- 763 )- 1 and g 2 = 0.560(1 + 21.1(fi 6 /i 2 ) 1 - 81 )- 1 . 
From the WMAP7 result (Komatsu et al. 2010), z ICC = 1091.3. The x 2 value of the shift parameter can 
be calculated by 

2 _ (i?_L725) 2 
0.018 2 

For the BAO observation from the SDSS spectroscopic sample of luminous red galaxy, the distance 
parameter is measured to be A = 0.469(n s /0.98)~ 0,35 ± 0.017 (Eisenstein et al. 2005), with the scalar 
spectral index n s = 0.963 from the WMAP7 result (Komatsu et al. 2010). The distance parameter can 
be expressed as 

2/3 



. l/2 -2/3 w i-l/3n-!/ 2 • o 1 / 2 f 

A=1] mo z bao £ ( z bao) sinn n k < / 

Jo 



dz 



(9) 



E(z 

where zbao = 0.35. The x 2 value of the distance parameter can be calculated by 

, _ (A- 0.467) 2 
Xbao - 0.017 2 ' 

The baryon mass fraction in clusters of galaxies from the x-ray gas (/ gas ) observation can be used to 
constrain cosmological parameters. On the assumption that the gas mass fraction in clusters is a constant 
and thus independent of redshift, Allen et al. (2004) obtained 26 observational / gas data. The baryon 
gas mass fraction / gas can be presented as 

- d SCDM (z)1 2/3 



/gas(z) — A 



(ID 



6a(z) 

where gJa = <^l/(1 + 2) 2 is the theoretical value of the angular diameter distance from cosmological 
models, d^ CDM is the angular diameter distance corresponding to the standard cold dark matter model 
(SCDM, fi M o = 1 for a flat universe), and A = [bVt h (2hf/ 2 ]/[(l + a)n M0 ], a = 0.19\/X, & is a 
bias factor motivated by gas dynamical simulations. The \ 2 value of cluster's baryon gas mass fraction 
(CBF) is 

2 [/ga?( 2 i) — /gas(Zi)] 2 
XCBF = 2^ 2 ■ ( 12 > 

i=l CT / g .= ^ 

The parameter A can be treated as a nuisance parameter by minimizing (Nesseris and Perivolaropoulos 
2007) 

X^bf = C- B 2 /A, (13) 
where A = EL/W/o7 gas ,*] 2 ' B = E \h**,if tf^v C = E [/gas,i/o-/ gM ,i] 2 , and / gaSji = 

[4 CDM (^)/^(^)] 2/3 . 

The Hubble parameter H(z) can be derived by 

1 dz 

H(z) = --— — . (14) 
1 + z dt 

From the Gemini Deep Deep Survey (GDDS; Abraham et al. 2004) observations of differential ages of 
passively evolving galaxies and other archival data (Nolan et al. 2003; Treu et al. 2001, 2002; Spinrad 
et al. 1997; Dunlop et al. 1996), Simon et al. (2005) have obtained the H{z) data at nine different 
redshifts (0.09 < z < 1.75). Recently, H(z) = 83.2 ± 2.1km/s/Mpc at z = 0.24, and H(z) = 
90.3 ± 2.5km/s/Mpc at z = 0.43 have been obtained by using the BAO peak position as a standard 
ruler in the radial direction (Gaztanaga et al. 2009). The x 2 value of the 1 1 H(z) data is 

2 N ^ [H ohs ( Zi ) - H{ Zi )] 2 
Xh= 2^ 12 • ( 15 ) 



=i 



'H,. 



The nuisance parameter H is also marginalized following the procedure used in calculating x 
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3 CONSTRAINTS FROM COMBINING GRBS, SNE IA, CMB, AND BAO 

In order to combine GRB data with the SNe la data to constrain cosmological models, we follow a 
simple way that avoids any correlation between the SNe la data and the GRB data (Liang, Wu and 
Zhang 2010): The 40 SNe points used in the interpolation procedure to calibrate GRBs are excluded 
from the Union2 SNe la sample used to calculate the joint constraints. Since the reduced 517 SNe la, 
42 GRBs, CMB, BAO, as well as CBF and H(z) are all effectively independent, we can combine the 
results by simply multiplying the likelihood functions. The best fit values for model parameters from 
the distance moduli of GRBs at high redshift obtained directly from the Union2 set, and SNe la, as well 
as the other joint observations (CMB+BAO+CBF+ iJ(z)) can be determined by minimizing 

2^2 2 2 ^2 ^2 

X = X M ,{42GRBs+517SNe} + XcMB + XbAO + XCBF + Xh ■ (16) 

In order to show the contribution of GRBs to the joint cosmological constraints, we also consider the 
X 2 value from the joint data (SNe + CMB + BAO + CBF + H(z)) without GRBs: x% = xf, { 55 7SNc} + 
Xcmb + Xbao + Xcbf + Xh> an d me j° mt constraints with GRBs + CMB + BAO + CBF + H(z) 
without the SNe la contribution is: x % = X^. {4 2GRBs} + Xcmb + Xbao + Xcbf + Xh- 

The joint confidence regions in {f2M<r^r c } plane with the combined observational data for the 
DGP model are showed in figure 1 . For comparison, fitting results from the joint data without GRBs 
are also given in figure 1. We present the best-fit value of {Omoj ^r c } with l-a uncertainties and the 
corresponding fix, as well as Xmin> Xmin/dof f° r the DGP model in Table 1. We also investigate the 
deceleration parameter for the DGP model. The deceleration parameter q(z) can be calculated by q = 
— 1 + (1 + z)E(z)~ 1 dE(z)/dz, where E(z) = H/Hq. And we could derive the transition redshift 
at which the universe of the DGP model switches from deceleration to acceleration (Zhu and Alcaniz 
2005; Guo et al. 2006) 

z T = -l + 2(-^) 1/3 (17) 

The best-fit values of zt of the DGP model are also summarized in Table 1 . 

With SNe la + GRBs + CMB + BAO + CBF + H(z), the best-fit values at l-a confidence level are 
{f^Mo, ^r c } = {0. 235lg 0.138^Q Q3g}; with the corresponding Qk ~ 0, which is near the line of a 
flat universe ((1 - n M0 ) 2 - 4ft rc = 0); with SNe la + CMB + BAO + CBF + H(z), the best-fit values are 
{fi M0 , O r J = {0.217±g;Jf|, 0.144i^35>; while with GRBs + CMB + BAO + CBF + H(z), the best-fit 

values are {f^Mo, ^r c } = {0-285lg 066j O-l^to 062}- Th ese results lead to more stringent constraints 
than previous results for constraint on DGP model with GRBs and/or other combined observations 
(Wang et al. 2009a; Wei 2010a; Bento et al. 2006; Davis et al. 2007; Reboucas 2008; Li et al. 2010). 
We also obtain the transition redshift zt = 0.6Tto '54 (ler) with the joint data including GRBs, which is 
more stringent and later the former result (zt = 0.86^0^) in Guo et al (2006). 

From comparing to the joint constraints with GRBs and without GRBs, we can see that the contribu- 
tion of GRBs to the joint cosmological constraints of the DGP model is a shift between the best fit values 
near the line which represents a flat universe, towards a higher matter density Universe (A^mo > 0). 
This situation has been also noted by Liang, Wu and Zhang (2010), and Liang, Wu and Zhu (2010), 
who comparing to the joint constraints with GRBs and without GRBs using the ACDM model, wCDM 
model, and Cardassian model. Also, a shift towards a later transition redshift can be found by comparing 
to the joint constraints of the DGP model with and without GRBs. It is shown that GRBs can give strong 
constraints on the DGP model when combined with CMB and BAO observations without SNe la, which 
has been also noted by Liang, Wu and Zhu (2010); Liang, Xu and Zhu (2010). 

4 CONCLUSION AND DISCUSSION 

In this paper, we investigate observational constraints on the DGP model including the cosmology- 
independent GRBs obtained directly from SNe la. Combining the GRBs at high redshift with the Union2 
set, the WMAP7 result, the BAO observation, the clusters' baryon mass fraction, and the observed 
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The DGP Model 




SNe+GRBs+Others 


SNe+Others 


GRBs+Others 


!T2mo 


u.zoo_ 074 
u.ioo_ 036 
0.033 


n 91 7+0. l^b 
u - zi ' -0.073 

U.144_ 035 

0.037 


U - ZOJ -0.066 
n i 99 +0.044 
U.1ZZ_ 062 

0.024 


Xmin 


595.95 


606.37 


77.93 


Xmin/d0f 


1.07 


1.09 


0.99 


ZT 


67 +u ' u:j 

u -°' -0.04 


n 74+o-ub 

u - '^-0.07 


U - Oi -0.16 



Table 1 The best-fit value of the DGP model parameters {riMo-^r c } an d with 1-cr uncer- 
tainties, Xm hl , Xmin/dof, as well as z T with SNe+GRBs+Others (CMB+BAO+ CBF+H(z j), 
SNe+Others, and GRBs+Others, respectively. 




0.2 0.3 0.4 0.5 



D.M0 

Fig. 1 The joint confidence regions in {fiMo-^r c } plane for the DGP model. The contours 
correspond to 1-cr and 2-cr confidence regions. The black solid lines, red dashed lines, and the 
blue dash-dotted lines represent the results of SNe+GRBs+Others(CMB+BAO+ CBF+if (z)), 
SNe+CMB+Others, and GRBs+Others, respectively. The black plus, red point, and blue star 
correspond the best-fit values of SNe+GRBs+Others, SNe+Others and GRBs+Others, respec- 
tively. The green line represents a flat universe which can be given by (1 — Omo) 2 — 4fi rc = 0. 



Hubble parameter data, we obtain {f^Mo, ^r c } = {0-235^0 ^74, 0.138^Q Q3g}, with the corresponding 
f^K = 0.033, which favors a flat universe. We also obtain the transition redshift of the DGP model 
zt = 0.671q'q4. These results breaks the degeneracies between the model parameters and leads to more 
stringent constraints than the previous results for constraint on DGP model with GRBs and/or other 
combined observations. It is shown that GRBs can give strong constraints on the DGP model when 
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combined with CMB and BAO observations. We conclude that GRBs could be used as an optional 
choice to set tighter constraints at high redshift on cosmological models. 

Zhu and Alcaniz (2005) tested the DGP model with the baryon mass fractions in clusters of galax- 
ies and the SNe la data to find that {Qmo^J = {0.29±g;gt, 0.21±g;g|}, and fl k = -0.36t o jg 
at 99.73% confidence level. Guo et al. (2006) also obtained a spatially closed DGP universe with 
fl k = -0.350^0 08° by using SNe la + BAO data. Zhu and Sereno (2008) used gravitational lens- 
ing statistics to find that the likelihood peaks at {f^Mo, ^r c } — {0.29, 0.12}, just slightly in the region 
of open models. These results seem to be in contradiction with the most recent WMAP results indicating 
a flat universe. However, constraints on the DGP model of the joint data including GRBs in this work 
are consistent with those obtained by Bento et al. (2006) using SNe la + CMB + BAO, and by Reboucas 
(2008) using SNe la + CMB; which favor a flat universe. 
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